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ABsIRACT: The Bruun Rule in shore crosion is presented systematically, starting from initial assump- 
tions and axioms towards a rigorous mathematical treatment. Mathematical treatment is given with detail 
to assist geologists in full understanding of the development. Cases of beach exeavation, littoral drift, ar- 
tificial beach nourishment and formation of a cuspate shoreline are treated by the theory, based on the 


Bruun Rule. 


INTRODUCTION 


The eoneept of shore erosion due to sea-level rise ex- 
pressed by Bruun (1962), later beeoming widely known 
as the “Bruun Rule” (Sehwartz 1967), belongs to the 
elass of eoneepts, whieh are enthusiastieally supported 
by some scientists, while immediately triggering a 
negative response from others. 

The ehronology of Bruun Rule studies (Fisher 1980a) 
ineludes laboratory and/or field observations, by 
Schwartz (1965, 1967), Dubois (1975, 1976, 1977a, 
1977b), Hands (1976, 1977, 1980), Fisher (1977a, 1977b, 
1977e, 1980b), and Rosen (1977, 1978a, 1978b, 1980); all 
of them strongly supporting Bruun’s econeepts. On the 
Other hand, the concept was eriticized by Swift (1976) as 
being of limited applicability. Another eritieism is that 
by Kaplin (1973), who is skeptieal about the validity of 
Bruun’s concepts; although supportively stating 
simultaneously that at least some of Bruun’s eoncepts 
were known by the Soviet school as early as 1946. 
Valuable comments on applicability of the eoneept are 
given by Gill (1979). 

lt appears that most of the controversy exists solely 
due to differences in interpretation of the Bruun Rule. 
The differenees arise, in their turn, due to lack of 
rigorous mathematieal formulation as well as lack of 
clarity in the statements of the initial assumptions of the 
Rule, Hence, our task in this paper is to present the 
Bruun Rule as a rigorous theory based on elearly stated 
assumptions. It is shown that the Bruun Rule has, in 
faet, a much wider field of application, than was 
Previously thought. The Bruun theory ean be applied 
for both “elosed” and “open” beach systems. We start 
from formulation of the Bruun Rule as in Schwartz 
(1965, 1967). 


THE BRUUN RULE AND ITS INITIAL ASSUMP- 
TIONS 


The Bruun Rule states the following (Fig. 1): 

l. A rise in sea level eauses erosion of the upper beach 
and shoreward displacement of the shore-water 
boundary. 

2. The change in sea level corresponds to translation of 
the transverse beach profile while retaining its 
original shape. 
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3. The material eroded from the upper beach is equal in 
volume to the material deposited on the nearshore 
bottom. 

4. The rise of the nearshore bottom is equal to the rise 
in sea level. 

5. The relationship between sea level rise a and 
shoreward displacement s of the beach profile is given 
by the formula (Bruun 1962): 


_ hs 
all 


(1) 

where: / isthe length of the transverse profile, 

h is the profile height, being the sum of sea 
depth at the distanee / from the shore and 
the shore elevation above the sea level. 


LIMITS OF THE BRUUN RULE 


As stated by Gill (1979) the rule only applies I, where 
there is sufficient energy; 2, where equilibrium has been 
attained; 3, where there is sufficient space in the subtidal 
area; and 4, if there is sufficient sediment. 


EXPERIMENTAL SUPPORT OF THE BRUUN 
RULE 


The two elasses of experiments, supporting the Bruun 
Rule, are laboratory (Sehwartz 1965, 1967) and field ex- 
periments. The advantage of laboratory experiments is 
that the “equilibrium profile’ whieh is praetieally 
unobservable on a beach, due to eontinuous variations 
in wave climate, ean be maintained in a laboratory 
where the wave elimate ean be set constant. 


LABORATORY EXPERIMENTS 


In the first experiment (Sehwartz 1965), the wave 
basin was 81.25 em wide and 115 em long with variable 
gradient of the bottom adjustable at 0°, 2.5°, and 4.5°. 
The waves generated had a period of 0.33 see +5%, an 
amplitude of 842 mm and a wave length of 15+ I em. 
Water depth ranged from 5 to 10 em, the sand used was 
a natural Ottawa sand, washed and sorted. 

It was found that 30 mins of wave attaek produced a 
beaeh profile whieh did not ehange with subsequent 
wave aetion, and this profile was assumed to be the 
“original profile” figuring in the Bruun statements. The 
water level was then raised by 10 mm, the wave 
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Fig. 1—Definition sketch of the variables in the Bruun Rule. 

The proto-profile AOB is translated towards the shore (positive 

x-direction) 10 its new position A’O'B’. However, no deposition 

occurs at the values of x<0. The segment DO is assumed 10 be 

the sharp boundary between the regions. In reality there will 

always be a smooth transition from the point D to the curve 
BO. However, this is neglected in the Bruun theory. 


generator was set up again and as soon as the 
equilibrium profile was reached water depth and depth 
of sediment at the outer edge of the shelf were 
measured. 

The profile retained its original shape after water level 
rise and was translated towards the shore, causing the 
shore recession. During the experiment the slope of the 
bottom was varied and it is remarkable that this change 
did not have any significant effect on water depth or 
sand heights accumulated at the lower portions of the 
beach. 

In the second laboratory experiment (Schwartz 1967) 
the wave tank was larger in size (100 x 232.5 cm) and the 
wave generator operated with variable periods. The sand 
was the same as in the previous cxperiment and water 
level was raised in different increments and at various 
rates. 

Profile shapes were observed before and after water 
level rise and the results indicated support for statements 
1, 2 and 3 of Bruun’s Rule. The most important state- 
ment (4), that the increase in sand height deposition on 
the outer shelf equals the rise of water level (hence, 
maintaining the constant water depth there), was sup- 
ported by the experimental results. 

These experimental results support the first 4 
statements of the Bruun Rule. However, formula (1) 
was not checked during the laboratory experiments. 
This was done in some of the following field ex- 
periments. 


FIELD EXPERIMENTS 


In 1964 an experimental program was established to 
study the variation of beach profiles under the variation 
of sea level due to neap and spring tides (Schwartz 1967, 
1979). At Herring Cove Beach in the Cape Cod National 
Seashore (USA) during a slight breeze, small waves ap- 


proaching the breaker zone were observed with 15-20 cm 
heights and periods of 3 sec. The direction of wave ap- 
proach during the summer is predominantly from the 
north-west. Shore drift at this time, is consequently, 
toward the south. Scdiment supply for this part of the 
Cape was provided by shore drift from the cliffs of 
glacial drift that form the outer coast. The nearshore 
bottom is characterized by a stcep slope prior to grading 
off to a gentler slope. 

At nearby Nauset Light Beach which bears 12° wes, 
of north and is exposed to the fetch of the Atlantic 
Ocean, small waves from the southeast approach the 
breaker zone on a calm day with a period of 8 secs and 4 
height of 30 cm. Erosion of the back beach cliff of 
glacial drift supplies scdiment that drifts northward. 
The nearshore bottom slopes gently to a bar 500 m 
offshore at low tide. 

The results of profiling during low and high water 
levels indicate parallel translation without significant 
change in shape during sea-level rise. The more dynamic 
regime at Nauset Light Beach means greater translation. 

Another field cxperiment, strongly supporting state- 
ment 2 of the Bruun Rule, was conducted by Dubois 
(1976) who measured nearshore profiles of Lake 
Michigan as water level seasonally rose from April to 
July, 1971. There were a total of 17 measurements of 
beach profiles during this period. The profiles presentyd 
by Dubois clearly indicate the parallel translation 
without significant changes in shape of profiles as water 
level rose. 

Using the Bruun formula (1), in which the value / was 
taken as a distance from foreshore base to the position 
of breaking waves, Dubois (1977a) calculated shore 
recessions and compared them with the observed values 
for various wave conditions recorded. The agreement 
between the values of shore recession observed in the. 
field and those calculated on the basis of Bruun’s for- 
mula (1) was remarkable. 

The above field observations were limited to only two 
sites and a few months in both Schwartz’s and Dubois’ 
experiments. The following field experiment by Hands 
(1976, 1980) was devoted to more long-term shore ero- 
sion. 

Hands (1976, 1980) described the response of Lake 
Michigan shores to increased water levcl over a 9 year 
period at 34 sites. He used formula (1) and demon- 
strated good agreement with observed erosion. His most 
interesting result was the clear transition zone between 
the area of bottom erosion and no erosion. Hands was 
the first to point out that actual length of the bottom 
profile (value / in eq. (1)) is not important, because it is 
the value A//, which is critical in the formula (1), and 
this value does not change much whether the length of a 
transition zone is considered or not. 

Observations by Fisher (1977a, 1977b, 1980b) of 
Rhode Island shoreline retreat over 35 years with a scope 
of 113 sites take into account an advance of sea-shore 
boundary due to a submergence (drowning) of the shore 
without erosion, separating this from shore retreat due 
to shore erosion, to which he applies Bruun’s formula 
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(1). The greatest contribution of Fisher is to consider the 
“open” system, in which not all the sediment volume 
eroded is deposited on the nearshore bottom. Sediment 
budget caleulations by Fisher require knowledge of the 
position of the point betwecn the beaeh erosion zone 
and the offshore deposition zonc (point Z on Fig. 1) to 
which various names had been preseribed (inflection, 
fulcrum, null point). 

It will be demonstrated below, however, that 
knowledge of position of sueh a point is not needed for 
cither “closed” or “open” systems, 

The longest projeet undertaken up to date was that of 
Rosen (1978b, 1980), who described the shore erosion 
On 146 beach units, due to sca level risc, over 100 years. 

The study area consisted of 350 km of estuarine 
Shoreline in the southern half of Chesapeake Bay. 
Bruun’s formula (1) was applied to caleulate an average 
Over 100 years’ shore recession rate and the results were 
compared to the measured values. 

On the basis of the above laboratory and field ex- 
periments it is clear that experimental support exists for 
the Bruun Rulc. We now havc to analyze the Bruun Rule 
Statements themselves. 


ANALYTICAL DESCRIPTION OF THE BRUUN 
EFFECT 

Statement 2 of the Bruun Rule, termed the “Bruun 
Effect” (Schwartz 1967), must be aeeompanied by state- 
ment 1, which indicates in what direction (onshore or 
offshore) the profile is translated during sca level rise. 
Considering these two statcments together as separated 
for the time being from the other Bruun statements, we 
Must introduce a notation @ for the rise of the profile, 
which is not required at this stage to be cqual to the sea 
level rise a. 

The bottom profile is considered as a function /(x,5) 
of the variable x and the parameter s. 

Assume that the initial shape of the profile, before 
erosion oecurred (s=0), was given by a eertain arbitrary 
function f,(x). Then the initial condition ean be written 
as 

Kx, s) = folx) (2) 

C=) 
The translated profile, shifted by s and lifted by a(s) is 
given by 

Sx, 5) = fx- 8) + uls) (3) 

Let us assume now, that we aeeept Bruun formula (1), 


equating @=a, which, when substituted into (3) gives the 
result: 


Axa) = fax- H +a (4) 


This is the analytical form of the Bruun effect, 
because, when a=0, the function f(x,a) reduces to the 
initial profile shape f(x); for any given positive sea level 
tise a, the profile is shifted shorewards in the positive 
x-direction by the value a//f and raised by the value a, 
retaining its original shape f(x), henee statements I and 
2 of Bruun’s theory follow from (4). (For the negative a 


(sea level fall), the profile moves offshore and 
downwards). 

Let us consider now the case in a ecrtain sense op- 
posite to rise in sea level, namely, elevation of the beach 
profile by the value b(s) due, for example, to artificial 
nourishment of a beach. Elevation of the beach profile 
while the absolute sea level remains still is equivalent to a 
fall of the relative sca level with respect to the nearshore 
bottom. Hence, according to Bruun’s statements 1 and 
2, the beaeh should retain its original shape, but must be 
translated seawards. 

If, once again, the Bruun formula (1) is utilized, with 
the value of profile rise b substituted instead of a, a for- 


mula analogous to (4) emerges: 


fx b) = fole+ PS +b © 


When therc is no beach nourishment (b=0), the 
beach profile reduces to the original profile shape /o(x). 

For any positive 6>0 (beach nourishment) the profile 
is shifted seawards by the value b//h, while retaining its 
original shape. (For b negative due, for example, to 
sinking of sediments offshore, littoral drift or bottom 
excavation, the profile is moving shorewards, i.e. ero- 
sion oeeurs.) 

If it happens that sea level rise a is aecompanied by 
beaeh nourishment, resulting in the bottom rise b, then, 
combining equations (4) and (5), the resultant beaeh 
transformation can be deseribed as: 


fla, b) = falx + + (b—a)) +(b+a) (6) 


where f(X) is the original shape of the beach. 

Note that the basie cquations (4), (5) and (6) were ob- 
tained solely on the basis of the first two statements of 
Bruun Rule, 


BRUUN RULE FOR CLOSED BEACH SYSTEMS 


Statement 3 of the Bruun Rule is the definition of the 
closed beach system, in which the volume of erosion is 
balanced by the volume of deposition and no exchange 
of beach material with the outer world exists. 

It may seem obvious that for the calculation of 
volumes of erosion and deposition one needs to know 
the position of the fulcral (null) point (see Fig. 1) and 
this need was, in fact, expressed by Dubois (1977a) and 
Fisher (1980). The difficulty here is that there eould be, 
in prineiple, several such fulcral points, as shown in Fig. 
2. Determination of their positions requires knowledge 
of beaeh profile shape and foreknowledge of the same 
values which one is going to ealculate on the basis of 
equating the volumes of erosion and accretion. 

This difficulty, however, can be by-passed if one 
prescribes opposite signs to the volumes of erosion and 
aeeretion according to the signs of differcnee in or- 
dinates between the protoprofile A(x) and the translated 
profile f(x,5). (Sign plus corresponds to aecretion; sign 
minus to erosion.) Then the definite integral taken along 
the x-axis will be equal to zero, when volumes (area be- 
tween the curves in Fig. 2) of erosion and aecretion 
balanee cach other. 
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Fig. 2—Multiple fulcrum points, corresponding to multiple 
areas of erosion (horizontal dashing) and deposition (vertical 
dashing). 


The important issue here is the limits of integration. 
As already mentioned, the value / (offshore beach 
length) is supposed to have been chosen in such a way 
that no erosion or deposition occurs offshore beyond 
this length, i.e. to the left of the point x=0 in Fig. 2. 
Therefore, it is reasonable to accept point x=0 as the 
lower limit of integration. On the other hand, if the 
beach profile retreat is equal to value s during sea level 
rise, it is reasonable to take the upper limit of intcgra- 
tion as (/+ 5s). Consequently, we conclude that Bruun’s 
statement 3 can be equivalently written analytically (in 
the co-ordinate system of Figs 1, 2) as 


t+ 


I= J Uf (xs) -f (9) dx=0 (7) 


o 

The first two statements of thc Bruun Rule were 
shown in the previous paragraph to Icad to the analytical 
expression (3) of the function /(x,s), in which the rise of 
bottom profile as a whole was denoted a(s). According 
to statement 4 we now equatc the rise of profile a(s) to 
the sea level rise a(s). Hencc, the function /(x,5), figuring 
in (3) can be written 


Ks) =o (x— 8) + as) (8) 


Let us emphasize that equations (7) and (8) comprise 
statements l-4 of the Bruun Rule, but the Bruun for- 
mula (1) (statement 5) is not used. On the contrary, we 
are going to demonstrate, that these four statements 
form the complete set of axioms, from which formula 
(1) can be derived. We start from several examples of 
profiles, having simple analytical expressions. 


ANALYTICAL PROFILES 
Parabolic profile 
Consider the profile, described by a parabola 
Jo(x) =x? (9) 


Let us note, that when x=/, the right end value of 
function (9) equals the profile height } (see Fig. 3), hence 
h= P. Substituting (9) into (8) and then into (7), one ob- 
tains 


Its Its 


f [(x +s)? + a(s) —x’]dx= Í [-2xs +$? +a(s)jdx = 


(10) 
— x35 +s?x+ xals)| = — Ps —ts?+(1+s)a(s) =0 
0 
Hence 
2 hs(1+ 5) 
as) = Esth vE (11) 
S Meg 


where the value was equatcd to A. Bruun’s formula (1) 
emerges rigorously from equations (7) and (8) for the 
parabolic profile. 

Now the position of the fulcrum point can be found 
as the point of intersection of the functions fo(x) and 
Axs) by solving the equation 
(notice that a(s) = hs/!=Ps/1= Is) 


x? =(x—s)?+ a(s) =(x—-5)? + ls 


O= —2xs +6 + ls 


s lia) t (12) 
2s 2 
The fulcrum point appears to be at the distance 
(/+5s)/2 from the coordinate system origin. 
Let us find the volume (area, as in Fig. 1) of acerction. 
This is the integral between the points x=0 and the 
fulcrum point x= “(/+s) 


Va (its) Va (its) 
VA =JS[(x-s)? +1s-—x7]dx=s{[(1+s) -2x]dx 
o o 
Ya (tts) 
=s[(/+5)x—x°] | ZASS) (13) 


als): -G($-) 
G= 13S 4-9 
a(s)= bs. Ga (4) 


=1+ Os + n(n-1)s° 
EREA 22 62 


EXPONENT 
f= e-t 


ARBITRARY WITH 
als) 
FLAT BOUNDARIES 
ARBITRA RYE > 
ANO SHOREWARD |a(s)= US+(&+a) 
LOPE £ (i+ 


Fig. 3— Analytical and non-analytical profiles with correspon- 
ding formula, relating sea level rise a to shore retreat s. 
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Analogously, the volume of erosion, i.e. between the 
fulcrum point x= '%(/+5s) and the shoreward limit 
x=/+5, is equal to 


VE= y [(x- s)? + is- x]dx= s{U+s)x—7}" 
Vaisa) Ya(tes) 


= —s(/+5)? (14) 


From this cxample it is seen that volumes of erosion 
and accretion are equal in absolute values, as it should 
be according to the Bruun statement 3. 

However, the actual volumes can be found ovi/y when 
the position of the fulcral point has been calculated. The 
impressive rcsult obtained was that the Bruun formula is 

` exact for thc parabolic profile. We shall see now, that 
this is not necessarily correct for other profile shapes. 


Cubic profile 
Let h= (15) 


then at the right end of the profile, where x =/, the value 
Jx) | TET: =/= h, the value A being the profile height in 
the Eat system of Fig. 3. Substituting (15) into (8) 
and then into (7), we arrive at 


= iG =s) + a(s) —x"]dx= -ws + = x75? — x5? + xa(s) [ 
=('s+ SPs! T La ts) a(s)=0 (16) 
Using the identity P= A, it follows: 
3s S nwad ys 
es ae 27 7 Seal) SAS SED S07 
K+ 5) 


One can recognize here the right-hand side of the 
Bruun formula, being multiplied by the correction Tac- 
tor G(s/f). When the ratio s// is small in comparison to 
unity, the value of the correction factor is also small. 
For example, when s//=0.1, the value G(s/)=1.05, 
Meaning that use of Bruun Tormula (1) for a cubic 
profile results in an error of 5%. For thc more realistic 
valuc s//=0.01, the factor G(s//) equals 1.005, meaning 
an error of 0.5%. We show below, that for other 
analytical profiles the errors could be much larger. 


Monomial profile 
EEO = xu (18) 
whcre n= tł. At the point x=/ the height of the profile 


is h=/, Substituting (18) into (8) and then (7), we can 
expand the integrand: 


n-i = 
~s) tals) -X= —nx s+” Oe S 


— WY 2) 34 2. (V's + a(s). (19) 


Upon integration, the powers of (/+s) appear, which 
arc cxpanded in the series such as 


U+ =r + Sy =r +n $+ MoD, Ft] 


-a S En 1) $ + 


n—1)(n-2) s? 
Aaa Nn 2), Bool 


Preserving only the terms with the value s// of power 
two and less, one obtains in the same manner as for the 
parabola and cubic function 


n @ a a & 
blt Z. $4 $y 


Bor K+ s/h 


= hs/l.G,„(s/) (20) 

lt is seen, that for large n the value of the correction 
factor can also be large and departure from the Bruun 
Tormula can be significant. For example, for n=10 and 
5/l=0.1 the valuc G, = 1.5. We shall discuss these results 
later and consider now another class of profiles, for 
which departure from the Bruun formula can also be 
large. 


Exponential profile 
Let fo(x) =e"—1 (21) 


The profile has a zero ordinate when x=0 and the 
height of the profile at x=/ is equal to h=e“'—-1, By 
varying the (positive) parametcr a various shapes can be 
obtained (Fig. 3). Substituting the translated function 


fas = et — 1 + a(s) = eve — 1 + a(s) (22) 


together with (21) into the integral (7) we obtain, with 
some simple manipulations 


1= | [e”(e™-— 1) + a(s) dx (23) 
= Lew =1) (ee=— 1) +(+ s)a(s) 
= 1 [ev —e*) + (1 —e°)] + (14 a9) =0 


As the value s is considcred to be small, we employ 
now the familiar Taylor expansions, retaining only the 
terms with s in the power two and less 


I-em=1-(l4+as ++ ...)2 - (as +25) 
I-ew=1-(1-as+%— ...)eas—% (24 
Substitution of (24) into (23) leads to 


=se" = 1) — Her + 1) + (14 sas) =0 


wherefrom, if one recalls that (e®—1)=h, the result 
emerges 


, _hs+ ale” + 1)s?/2 3 
S (25) 


It is scen that for œ small in comparison with unity the 
Bruun formula follows once again. However, for large œ 
(sce Fig. 3) the second term in the numerator may 
become predominant and it grows alarmingly fast with 
the increase of a. We again temporarily postpone the 
discussion of the reasons for such large errors and in- 
stead return back to the case of the linear profile, which 
we avoided discussing before, quite deliberately. 
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Linear profile 


Let fo(x) = kx; flx,s) = k(x— s) + als) (26) 
Where the slope k = a 
The integral (7) then reduces to: 
= Tike- s) + a(s) — kx]dx= aint —ks+a(s)\dx 
=- ks +a(s)=0 š (27) 
Hence, 
a(s) =ks=lis/l (28) 


and the Bruun formula emerges, exactly. Notice, that in 
this particular case of lincar profile (and in this case 
only) the integrand in (27) equals zcro identically, in- 
dependently of the value of x, as soon as (28) is ued. If 
we recall the case of the parabolic profile, where to find 
the fulcrum point we need to equalize the integrand to 
zero, it becomes obvious that for a linear profile every 
point xv is the fulcrum point; hence no crosion or accre- 
tion really occurs. 

As any analytical function f(x) can be expanded as a 
power serics, it follows from equations (11), (17), (20) 
and (28), that Bruun formula (I) is rigorous only for 
generalised parabolic profile 

SUX) = CX + C2X? (29) 
where the constants ¢, and c, are arbitrary. For higher 
order of profiles the Bruun formula (1) holds only as an 
approximation, as summarised in Fig. 3. 


Non-ANALYTICAL PROFILES AND BOUNDARY CONDITIONS 


One may notice that when the profile fox) is shifted 
shorewards by the value s, the portion of the prokile 
which originally corresponded to the values x<0, ap- 
pears at the positive part of the x-axis (see Figs I, 2). 

This means, that it is insufficient to know the profile 
shape fo(x) strictly within the interval O<x=s/, to 
calculate the integral (7) one nceds also to know the 
behavior of the function fo (x) beyond the points 0; /, 
namely, in the wider interval { —s;/+ s]. 

Hence, the outer boundarics of the segment {0;/] are 
important, and the shape of the profile at these boun- 
daries —ssx<0 and /<x</+s forms some kind of 
boundary condition (not to be confused with the boun- 
dary conditions for differential cquations). 

We are going to demonstrate now that these arc 
precisely the boundary conditions, which determine 
validity or invalidity of the Bruun formula (1). 

Let us consider an arbitrary function Ax) with its do- 
main of definition in the interval [-s<x</+s]. (The 
subscript zero is dropped out here for simpncity of nota- 
tion only.) 

Using the general formulae (7) and (8), th: integral 1 is 
as follows: 


1= Jifx—s) + a(s)—fWjdx = lpse lla =l (ea) 


where 


ts 


l= ‘j fx-s)dx; h= ‘Vals)dx; l; = [Avid (30) 


0 
By a change of variables the integrals may be rewritten 
as: 


L = |x- s)dx= {faydx= { Siodx+ {fdr= (31a) 


Siodx+ ÍJ (-—x)dx. 


~ oe © 


L= | a(s)dx=(I+ s)a(s) (31b) 
L= | fde= i fydx+ ffoydx= | Sixydx+ 
[fix+ Ddx Glo) 


Combining the expressions (31a, b and c) according ta 


t 
(29) it is seen, that the integrals f(x)dx cancel each other 
and: i 


s 


1=|{((-/)—fwt Ddx+ (+ sals) (32) 


According to cquation (7), the value 1 should equal 
zero and the rigorous expression for the relationship a(s) 
follows: 


aij = Sanna le (33) 


i+s 


Returning back to the temporarily dropped notation 
fax), it is convenient to rewrite (33) in the following 
form: 


a(s)= ne F(s)= Sox + N=-fol —x)\dx 

The main advantage of (34) is that it makes clear the 
role of boundary conditions. It is scen from (34) that - 
behavior of the profile /,(x) outside the interval {0; 4] 
contributes to the value F(s), while shape of the profile 
within this interval is of no importance at all. 

To clarify this statement, consider the two classcs of 
non-analytical profiles. 

(1) Consider onc class of profiles, which are arbitrary 
between the points [0, A and horizontal beyond these 
points (Fig. 3). The profile is assumed to have the height 
h, then 

Slx+ N=fol) =h 

Sol — x) = fo(0) = 0 (35) 

Calculation of the function F(s) in (34) gives, after 
substitution of (35) in it 


(34) 


F(s)= | hdx=hs (36) 
Hence, for such a profile 
hs hs 

or O N 7 

tae +p ey 


and Bruun formula (1) follows, if the value s//is small in 
comparison to unity. 

(2) Consider another class of profiles, which again are 
arbitrary between the points [0,/], but have extensions 
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beyond the terminal points as linear funetions (Fig, 3) 
Jax +N = h + Bx—shoreward 


fol — x) = -axy -— seaward (38) 

The parameter a is, in fact, the seaward slope of the 
profile, while 8 is the shoreward slope. 

Substitution of (38) into (34) gives 


` 


-2 
F(s) = | (h+ 6x+ axjdx=[hs+ (84 IE (39) 
0 

It is seen, that the seaward slope œ is aeeompanied in 
formula (39) by the corresponding shoreward slope 8 
and both slopes play equal roles. 

Partieularly, when the slopes are horizontal, (i.e. 
a= 8=0) formula (39) reduces 10 (36). 

Another case of reduetion to (36) occurs when 
a=—, i.e. seaward and shoreward slopes are of op- 
posite signs. 

An interesting case emerges when slopes a=6=k, 
where k=H//. Then it follows from (39) 


F(s) = hs + ks? = hs(1+ 7) (40) 
and the Bruun formula appears rigorously once again 
S S 
hs(t+ 7) _listt+ = ist 


as) = 
(as Kd Ge 


(41) 


Henee, when both seaward and shoreward slopes of 
Otherwise arbitrary profiles are equal to the value 4// of 
the profile, the Bruun formula (1) is exaet. 

The results obtained so far for the profiles of various 
shapes are depicted on Fig. 3 with the corresponding ex- 
act formula a(s). 

One can see from Fig. 3 that steepness of monomial 
and exponential profiles grows with the inerease of their 
Order and so does the error in the use of the Bruun for- 
mula (1), Hence, the boundary conditions, rather than 
the shape of the profiles, are responsible for accuracy of 
the Bruun formula, It is seen from (39), that the value 
F(s) ean depart from the Bruun value fs if the seaward 
and/or shoreward slopes are large. Reversely, one can 
always estimate the error, given by the Bruun formula, 
by means of exaet formula (34). 

Consequently, we may eonclude that Bruun’s formula 
(1) follows from the exaet formula (34) as a very ae- 
curate approximation, providing the seaward and 
shoreward slopes of the profile are not too steep. In 
turn, formula (34) is based rigorously on equations (7) 
and (8), which are the mathematical formulation of the 
first four statements of the Bruun rule. 

We show that the same Bruun statements form the 
basis for the theory of erosion/aceretion in open beaeh 
systems. 


BRUUN RULE FOR OPEN BEACH SYSTEMS 


It is convenient to define an open beaeh system as one 
in which an additional volume (+ V) of sediment is sup- 
plied from outside the system (for example, beaeh 
nourishment and/or littoral drift) or, alternatively, 
some volume of sediment (— V) is removed from the 


system (for example, due to exeavation of the nearshore 
bottom and/or to littoral drift). 


BEAC!) NOURISHMENT 


Consider the ease of a positive volume (+ V) of sedi- 
ment being supplied to the system, while the sea level a is 
still. This case was partly touched in section 4 where 
equations (3) and (4) were derived solely on the basis of 
the first two of Bruun’s statements. 

In this section we denote the horizontal profile 
displacement due to supply or removal of sediment as r, 
leaving the notation s for horizontal profile 1e1reat solely 
due to sea-level rise. Correspondingly, equation (3) in 
the new notation is 


Ax) =folx +r) + bir) (42) 


where D(r) is positive, when rise of the profile occurs due 
to income of positive volume (+ V) of sediment into a 
beaeh system and horizontal displaeement r due to the 
same factor is taken as positive for seaward displaee- 
ment of the profile. 

The increment in volume of sediment in a beaeh 
system due to translation of the profile given by (42) can 
be ealculated as an integral 1, analogous to (29). 

However, the limits of integration must be different 
now, As the profile is shifting seawards (positive r), part 
of the shifted profile will appear /eft of point x=0 (Fig. 
4), beyond the limiting length / offshore. But our initial 
assumption was that a profile is undisturbed to the left 
of point x=0. Henee, the lower limit of integration 
should be taken as x=0. 

On the other hand, as the profile is shifted seaward, 
no disturbance of the profile oecurs beyond point x=/, 
which shoufd be taken now as the upper limit of integra- 
tion. Consequently, the additional volume (+ V) of sedi- 
ment supplied to the beach should be equal to the in- 
tegral (we again temporarily drop the subscript zero in 


SAX) 
l=] [fv+)+b()-/W]dx=1,4+1,-L=V (43) 


o 
where 


Lh= § fixtndx= ee ee d (44a) 
L= | b(ndx= 1b) (44b) 
L= { fiodx. (44c) 


Adding to 1, eq. (44a) the self-eancelling pair of in- 
tegrals 


F fdx- dx (45) 
0 0 
we ean rewrite l, as 


R= f fx)dx+ Í fodx+'{foodx— Í fx)dx (46) 


where the sum of the first two integrals exactly equals lz, 
thus eaneelling it. 
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STILL SEA LEVEL 


Fig. 4— Deposilion of sediments, resulting in rise of profile by 
the value b(r) and seaward shift by the value r (r is positive 


seaward). 
Consequently 
ler r 
l= {[fiddx—- | fOddx+ lb(n=V (47) 
t 0 


and the formula for b(r) emerges (returning back to 
subscript zero) 

bin = 4(V— Vfotssdx+ | fod) (48) 

One has to notice that the profile shape between the 
points x=rand x=/ does not figure in (48) and value of 
the integrals is determined by the behavior of the profile 
shape in the neighborhood of its terminal points .=0 
and x=/. Let us assume that the profile is nearly flat 
around x=0 and x=/, i.c. at the interval [0,7] the func- 
tion (xX) =0 and at the shoreward end f(x) = / at the in- 
terval [/, /+7r], then the first integral in (48) reduces to 
the value Ar and the second to zero. Formula (49) then 
emerges as: 


b(n) =(V=hr)/1 (49) 


It can be shown, by analogy with the previous section, 
that (49) is a very good approximation of the exact result 
of (48) when the boundary conditions are such that the 
seaward and shoreward slopes are not too steep and the 
difference in their slope angles is small. 

The expression for r follows from (49) 


r=(V—bD/h - (50) 


On the other hand, we did have another expression 
for the same value r in equation (5) 


_ bl 


7 (51) 
Combining (50) and (51), we obtain 
ae ee 
r= 5, A b= 5, (52) 


Consequently, equation (42), becomes dependent on 
the volume of sediment supplied to the beach 


fix Y=flx+ 5 (53) 


We want to emphasize, that the values r and b in (52) 
and formula (53) are independent of the profile shape. 


If the beach loses sediment (due to cxcavation of the 
bottom, littoral drift, etc.), then the value V in (52) and 
(53) should be taken as negative and the beach profile 
shifts shorewards and downwards, because the values / 
and b become negative. 

When beach nourishment (excavation) is accom- 
panied by sea level rise, it follows from equations (6) 
and (53), that the new beach profile is described by the 
equation 


fia Wafixt LE -at de tah, 54) 


where positive valucs V and a correspond to beach 
nourishment and rise in sca level. 

Equation (54) was obtained, after all, solely on the 
basis of the first four Bruun Rule statements. 

The previous treatment was performed for a two- 
dimensional cross-section ol a beach, which may also be 
considered a slice (of unit thickness) of an actually three- 
dimensional beach. 


BRUUN RULE FOR THREE-DIMENSIONAL BEACH SYSTEM 
WITH LITTORAL DRIFT 


Consider the three-dimensional portion of shoreline 
in the coordinate system of Fig. 5 (note that the y coor- 
dinate is now directed alongshore). Any vertical plane, 
parallel to the plane ZOX, dissecting the shore, forms in 
its cross-section a two-dimensional profile, such as 
described in the previous sections. However, there is no 
need to assume that all the cross-sections have the same 
or similar shape. For each cross-section with a longshore 
coordinate y equation (54) holds; therelore independent 
of the shape of the particular cross-section, its seaward 
motion (denoted R(y) from now on) and vertical uplift 
(denoted @()) is equal to 


RO= nos YW) — ayy] (55a) 


ay) =a + [VO)2iw)) (55b) 


where a is the rise in sea level as before and the functions 
hO); 10); VO) are dependent on the longshore coor- 
dinate y. Consider the equation of continuity or the 
volume conservation equation for longshore sediment 
drift (see, for example, Komar, 1976) 


aV/at= —aQ/ay 


where Q is the longshore volumetric rate of sediment 
transport (measured, for example, in m?/day). 

The physical sense of equation (56) is quite obvious: it 
states that the time-rate (speed) of volume of sediment 
removal (sign negative in 56) from a beach equals the 
longshore rate in change of volumetric transport along 
the beach. This is illustrated in Fig. 5, where both func- 
tion O(y,f) and dQ/dy are plotted along the longshore 
coordinate y for one specilic type of dependancy Q(y). 

It is reasonable to assume further that the basic 
parameters /i(y) and /(), figuring in (55), do not change 
much in time due to shore erosion (accretion), and only 
the values, R, 4, a and V in (55) are variable in time. 
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Fig. 5— Definition sketch for longshore littoral drift. Axis y 
directed alongshore. Profile shapes could be different for 
various y, such as y, and y. 


Consequently, equations (55) can be rewritten to incor- 
Porate the above assumption as 


— Laren 
EO) =o) tet) a(t)i(y)] (57a) 
Aly,t) = a(t) + yee (57b) 


Notice that sca level rise a in (57) is assumed to be 
dependent on time, but independent of longshore coor- 
dinate y, meaning that extremely long coastlines, where 
sea-level varies along the coastline, are not considered. 
Equations (57) give the complete steady-state solution of 
the problem of beach nourishment (excavation) in the 
Presence of variations in the time sea level rise a(r) for 
any distribution of V(y, along the shore and in timc, 
providing the values A(y) and /(y) are known. 

Differentiating (57) with respect to time ¢ one obtains 


OR/at= Y2dV/at —I)da/dt/hO) (58a) 

ðā/ðt = ða/ðt + l(y):dV/at (58b) 

Substituting instead of dv/dr in (58) the equal value 
~—dQO/ay from (56), it follows 

ðR/ðt = —[%2°dQ/dy+ Ky):da/dg/hiy) (59a) 

da/at = da/at — YI(y):dO/ay (59b) 


Integrating equations (59) with respect to time, one 
obtains the final formulae for calculation of shoreward 
displacement R and profile risc @ as functions of 
longshore sediment transport Q and variations in the 
time of sea level rise a(r) 


ROD = = plg [2Par+ 10). ao (60 
TY) = a) -ay | ar (60b) 


We illustrate the use of (60) by two examples. 


Example | 


Let a()=0 and Q0,)=y sin wf, meaning that the 
sediment transport rate is linearly distributed along the 
shoreline and varies periodically in time. 

Then 


0007, 0/dy =sin wt (61) 


and from (60) one obtains by integration: 


ROTS -as wt +C 


-agg COS wt + D 


(62a) 


ay, = (62b) 


where C and D are the constants to be determined from 
initial conditions. Assuming them as R(y,)=0 and 
ay, í) =0 (no initial crosion—accretion) when 4=0, it is 
easy to find from (62) 


jl ee l 


hence, the seaward shore displacement and rise of the 
bottom are: 


RVD =at ~cos wl) 


aty,) = slay" — cos wl), (64b) 


meaning periodic (for example; seasonal) variation in 
shore erosion, dependent on the original longshore 
shape (functions /A(y) and /(y)) of the coastline. 


(64a) 


Example 2 


Let a(t)=0 and a tropical cyclone (hurricane) moves 
alongshore, causing a propagating wave-like longshore 
transport O(y,/ with magnitude Qs 


OU, À) = Qo cos (Ky + wt) (65) 
Then 
ga= — KQ» sin (Ky + wl) (66) 


and upon substitution of (66) into (60a, b) and integra- 
tion 


RO,D= — ee cos (Ky + wf) + CO) 


RA (67a) 
TERHI- = Re cos EuN + DO) (67b) 


2wl0) 


Here instead of constants C and D, as in the previous 
example, the functions CO) and DO) should be found 
from initial conditions. One obtains from (67) for t=0 


a KQo, o 
CU) = Toho) =——cos Ky (68a) 
Dy) =S Ky (68b) 


and the final result emerges from (67) upon substitution 
of (68) 
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RO,O= MO [cos Ky — cos (Ky + wt)] 


2wh(y) (ey 
ay, n=K Qetcos Ky—cos (Ky + wt) (69b) 


2uwl(y) 


Note the appearance of the time independent term cos 
Ky in (69), meaning that after passage of a cyclone (hur- 
ricane) an irreversible cuspate change of a coastline oc- 
curs, if the rate of sediment transport could be approx- 
imated as a propagating wave (65). 


DISCUSSION 


The first four statements of the Bruun Rule form a 
non-contradictory set of axioms, which need no other 
axioms for development of the theory of shore erosion; 
incorporating not only sea level rise, for which the 
statements were originally intended, but also the prob- 
lems of beach nourishment, excavation and littoral 
drift. 

The most important of all assumptions of the theory, 
its cornerstone, is the statement, that beach profile 
essentially preserves its shape during process of shore 
evolution, be it due to sea-level rise, wave attack, beach 
nourishment etc. Of course, it is recognised, that some 
temporal variations, perhaps caused by a sudden storm, 
can be superimposed on the conservative shape of beach 
profile, but these variations are usually short-lived and 
can be discarded in the first approximation. 

The experimental data available support this assump- 
tion. No doubt, further experimental work is necessary. 

The most important limitation of the above theory is 
its static nature. 

This means, that the beach system is considcred only 
in two states: initial state, before evolution started and 
final state, after evolution was completed. The transition 
process between these two states is beyond the scope of 
the present theory. Clearly, this shortcoming must be 
rectified. This is our direction for future research. 
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